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ABSTRACT 

In  order  to  relate  igniter  energy  output  to  motor  ignition  a  model 
of  solid  propellant  ignition  is  given.  Difficulties  inherent  in  the  non¬ 
linear  partial  differential  equations  of  heat  and  mass  diffusion  describing 
the  ignition  system  prevent  the  model  from  being  completely  analytical. 
However,  certain  results  are  derived  from  the  equations,  and  when  combined 
with  arc-image  furnace  ignitabillty  data,  a  reasonably  coherent  view  of  the 
process  is  obtained.  Numerous  examples  of  ignitabillty  data  are  given. 

Analysis  of  heat  and  mass  transfer  from  igniter  to  motor  grain  is 
given  for  the  case  of  propellant  type  igniter.  This  analysis,  along  with 
the  pertinent  ignitability  data,  is  applied  to  Ignition  of  the  2nd  stage 


INlnuteman  motor. 


I. 


INTRODUCTION 


In  order  to  predict  performance  and  specify  reliability  of  solid 
propellant  rocket  ignition  systems,  it  has  been  found  necessary  to  have  n 
model  of  propellant  ignition.  In  other  words,  rational  design  of  an  ig¬ 
nition  system  implies  that  igniter  output  be  based  on  ignition  energy  and 
pressure  requirements  for  a  particular  propellant,  perhaps  having  specified 
surface  preparation,  or  having  undergone  prescribed  environmental  con¬ 
ditioning. 

Attempts  to  establish  ignition  energy  requirements  have  led  to 
various  models  of  the  ignition  process.  Until  recently,  a  thermal  model 
of  ignition,  including  heat  generation,  yielded  the  usual  methou  of 
analysis.  This  model,  which  is  adequately  described  in  References  (l) 
through  (5),  is  essentially  based  on  the  assumption  of  solid  phase  reactions 
usually  taken  as  first  order.  Gas  phase  reactions  are  neglected.  Using  this 
model,  investigators  were  able  to  explain  data  acquired  from  ignition  tests 
where  hot  wire,  hot  plate,  and  hot  bath  apparatus  wore  used  as  heat  .'•■ources. 
These  types  of  apparatus  virtually  exclude  the  ambient  atmosphere,  a.nd  so 
diffusion,  as  an  ignition  parameter. 

In  the  ignition  analysis  given  here,  the  thermal  model  of  ignition 
is  taken  as  the  starting  point.  Results  found  by  investigators  applying 
this  model  will  be  used  liberally.  An  extension  of  the  thermal  model  will 
be  made  in  order  to  account  for  known  pressure  sensitivity  when  ignition 
occurs  in  a  gaseous  medium.  The  problem  of  heat  and  mass  diffusion  in  a 
reacting  gas  is  found  to  be  extremely  complex,  but  an  approximation  can  oe 
made  which  will  yield  some  information. 
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Test  data  obtained  by  means  of  the  Stanford  Research  Institute'  arc- 
image  furnace  are  given.  This  ignitability  data  has  been  obtained  for  a 
variety  of  propellant  compositions,  having  several  surface  conditions. 

Finally,  igniter  energy  and  pressure  output  is  defined  by  means  of 
mass  and  energy  balances  applied  to  a  propellant  ignition  system.  The 
entire  analysis  is  applied  to  the  ignition  of  the  Mlnuteman  2nd  stage 
motor.  An  ignition  safety  factor  is  defined  and  computed  for  the  stage 
in  the  ignition  train  consisting  of  motor  ignition  from  main  igniter  charge. 
II,  PROPELLANT  IGNITION 

Recent  ignition  research  has  demonstrated  that  the  ignition  process 
is  pressure  dependent  (Reference  (6),  (7)  and  (8)),  Hence  any  ignition 
model  must  consider  both  propellant  surface  phenomena  and  vapor  phase 
reactions  in  the  vicinity  of  the  propellant  surface.  In  the  model  presented 
here  it  will  be  assumed  first  that  the  solid  propellant  undergoes  a  primary 
decomposition  from  solid  into  gases,  which  undergo  further  reactions  before 
the  final  combustion  products  are  reached.  The  primary  decomposition  is 
a  result  of  a  heat  flux  at  the  propellant  surface  which  is  sufficiently 
high  so  that  o.nly  a  very  thin  layer  of  propellant  is  affected  before  ig¬ 
nition  occurs.  In  this  case  heat  generation  by  solid  phase  reaction  need 
not  be  considerea.  Secondly,  gas  phase  reactions  of  decomposition  products 
are  assumed  to  occur  which  will,  under  good  ignition  conditions,  transfer 
a  sufficient  amount  of  energy  back  to  the  propellant  surface  to  sustain 
the  combustion  of  the  solid  propellant. 

The  time  interval  oetween  ig.niter  flux  onset,  when  the  propellant 
surface  is  first  exposed  to  the  igniter  heat  flux,  and  the  instant  when 
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surface  decomposition  begins,  is  called  the  thermal  induction  Interval.  It 


is  characterized  by  accumulation  of  heat  in  the  propellant,  establishment 
of  a  temperature  gradient  at  the  exposed  surface,  and  finally  surface  de¬ 
composition. 

The  time  interval  between  commencement  of  propellant  surface  decom¬ 
position,  under  the  action  of  the  igniter  induced  heat  flux,  and  self- 
sustaining  combustion,  is  called  the  chemical  Induction  interval.  It  is 
characterized  ay  an  accumulation  of  active,  gas  phase  species  in  the  vicinity 
of  the  propellant  surface,  further  exothermic  reactions  of  the  species,  and 
finally,  under  good  ignition  conditions,  a  reaction  rate  sufficiently  higii 
to  cause  self-sustained  combustion  of  the  solid  propellant.  It  must  be 
emphasized  here  that  time-to-ignition  will  be  the  sum  of  thermal  and  chemical 
induction  intesTvals, 

The  two  pertinent  ignition  time  intervals,  thermal  and  chemical 
induction,  will  be  treated  separately.  It  will  be  found  that  the  value  of 
the  thermal  induction  interval  can  be  calculated  with  resonable  certainty. 

In  contrast,  the  only  information  gained  on  the  chemical  induction  interval 
will  be  its  dependence  on  such  parameters  as  ambient  pressure.  The  reasons 
for  this  are  the  complexity,  and  non-linearity,  of  the  partial  differential 
equations  governing  mass  and  energy  transfer  in  a  gas  undergoing  exothermic 
chemical  reactions.  Analytic  solution  of  tne  equations  appears  impossible, 
and  numerical  solution  usually  obscures  significant  trends.  Fortunately, 
data  made  available  through  arc-image  furnace  ignition  tests  determine  time 
to  ignition  at  a  variety  of  ambient  pressures.  The  thermal  induction  in¬ 
terval  may  be  calculated  and  subtracted  from  the  test  data,  with  the 
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difference  being  the  chemical  induction  interval.  Combining  this  experi¬ 
mental  data  with  the  analysis  in  this  section  permits  extension  of  the 
data  so  that  the  chemical  induction  Interval  can  be  obtained  for  other 
significant  ignition  parameters.  That  is  to  say,  combining  the  meager 
analytical  results  with  experimental  data  from  ignition  tests  will  result 
in  a  fairly  comprehensive  picture  of  gas  phase  ignition. 

A,  THERMAL  INDUCTION  INTERVAL 

'Vhen  solid  rocket  propellant  is  exposed  to  a  high  intensity 
heat  flux,  as  in  the  case  of  normal  ignition,  the  surface  of  the  propel¬ 
lant  is  assumed  to  undergo  a  primary  decomposition  from  solid  into  gases. 
These  primary  decomposition  products  are  assumed  to  react  further  before 
final  decomposition  products  are  reacned.  In  normal  ignition,  the  igniter 
induced  heat  flux  is  sufficiently  high  so  that  only  a  very  thin  layer  of 
propellant  is  affected  before  ignition  occurs  (References  (1)  and  (2)) 
and  in  this  case  heat  generation  by  solid  phase  reactions  need  not  be  con¬ 
sidered.  The  time  required  to  bring  about  primary  surface  decomposition, 
which  is  identical  to  the  thermal  induction  interval .  is  calculated  by  a 
simple  thermal  analysis  at  the  propellant  surface.  Use  of  a  single,  con- 
start  ignition  temperature  is  based  on  results  found  by  Altman  and  Grant 
(Reference  (3))" 

As  an  illustration  consider  the  case  of  an  igniter  with  a 
nearly  constant  mass  flow  rate,  and  hence  a  nearly  constant  heat  flux. 

When  the  gas  temperature  of  the  igniter  combustion  products  is  high  the 
heat  flux  at  the  propellant  surface  is  nearly  independent  of  surface  tempera- 
tur.’.  In  this  case  the  transient  temperatures  inside  the  propellant  grain 
are  given  by, 
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where 


thickness , 


T(x,t)  - 


The  dimensionless  numbers  are, 


X  =  x/Z^fat 

T(x,t)  =  temperature  at  position  x,  and  time  t 
=  initial  ambient  temperature 
q  =  igniter  induced  heat  flux 
a  =  propellant  diffusivity 
k  =  propellant  conductivity 

Now  let  t^  be  the  time  required  to  raise  the  critical  layer 
x^,  to  the  ignition  temperature,  Then  from  Eq.  (l) 


Equation  (2)  then,  establishes  relations  between  the  thermal 


induction  interval,  t^,  the  ambient  temperature,  T^,  and  the  various  pn'- 
pellant  and  igniter  properties. 


B.  CHEMICAL  INDUCTION  INTERVAL 

The  chemical  reaction  rate  will  be  taken  as  first  order  which 
is  in  agreement  with  a  burning  rate  exponent  of  one-half.  The  equations 
of  mass  and  energy  conservation,  for  an  immobile  system,  become  respectively, 


where 


d  C  _  ^  ^i^C  „  jr -E/RT 
d  t  -  ^  ^  -  ^1^" 


d  T  d^T 


-E/RT 


C  =  concentration  of  reactants 


(3) 

(4) 


=  pre-exponential  factor 
E  =  activation  energy 
R  =  universal  gas  constant 
T^  =  explosion  temperature 
D  =  mass  diffusivity 
tC  =  gas  thermal  diffusivity 
p  =  mixed  gas  density 
X  =  distance  from  propellant  surface 
The  concentration  of  reactants,  C,  is  defined  in  terras  of  its 
partial  pressure,  by 


C  =  (Pj^/bT) 


where  b  is  the  specific  gas  constant. 


In  order  to  obtain  a  rough  approxination  to  the  reactant  weight 
density  buildup  before  ignition  the  mass  balance,  Eq.  (3),  will  be  solved 
assuming  constant  temperature.  In  this  case,  Eq.  (3)  and  the  boundary 
conditions  are, 


62^ 
d  2  ~ 

X 


(5) 


=  Z^e 


-E/RT 


,  constant  mass  flux 


q»  fay^ 


C=0,  x2  0,  t=0 

The  mass  fluxj'p,  is  approximately  related  to  the  heat  flux, 


where  c  is  the  heat  capacity  of  the  solid  propellant,  r  is  the  linear  re¬ 
gression  rate  of  the  propellant  surface,  and  the  propellant  density. 
Existence  of  such  a  mass  flux  prior  to  ignition  has  been  verified  by  arc- 
image  furnace  ignition  tests  (Reference  (6)). 

The  solution  to  Eq.  (5)  at  the  surface,  x  =  0,  is. 


C  =  7-^  erfyi^  (6) 

V  K^D 

Before  ignition  occurs,  the  temperature  rise  near  the  surface 


should  be  small.  Then  the  reaction  rate  varies  directly  with  concentration, 
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Using  this  relation,  the  pressure  dependence  of  the  chemical  Induction 
interval  can  be  established,  and  is 


where  W  is  defined  by  Eq.  (7). 

The  plot  of  I  J  against  the  inverse  square  of  the  right 
hand  side  of  Eq,  (7),  which  is  shown  in  Figure  (1),  shows  that  for  a 
sufficiently  low  pressure  the  chemical  induction  interval  approaches  infinity. 
In  other  words,  for  a  given  propellant  there  is  a  pressure  below  which  ig¬ 
nition  will  not  occur.  Furthermore,  for  moderate  pressures  the  chemical 
induction  interval  approaches  zero.  Inspection  of  test  data  shows  that  for 
ambient  pressure  exceeding  about  six  atmospheres,  the  chemical  induction 
interval  is  nearly  non-existent  for  polyurethane  composite  propellant.  By 
conducting  ignitability  tests  at  high  pressure,  values  of  the  thermal  in¬ 
duction  interval  are  isolated,  and  by  using  steady  state  values  of  the  ig¬ 
nition  temperature,  the  value  of  the  critical  layer  thickness  is  obtained. 


o  T 


Ignitability  tests  conducted  near  the  lower  ignition  pressure  threshold 
level  will  allow  determination  of  the  pressure  parameter,  VT,  (Eq.  (?)). 

An  additional  ignitability  determination  in  the  pressure  range  betw'een  zero 
and  infinite  chemical  induction  interval,  will  suffice  to  fix  the  entire 
ignitability  curve  with  some  degree  of  certainty. 

C.  IGNITABILITY  DATA 

Ignitability  measurements  have  been  made  on  a  number  of  pro¬ 
pellant  compositions,  under  several  conditions  of  heat  flux  and  environ¬ 
mental  pressure.  The  device  used  to  acquire  the  data  was  an  arc-image 
furnace,  operated  by  the  Stanford  Research  Institute.  Propellant  samples, 
in  a  fixed  pressure  environment  of  gaseous  nitrogen,  are  tested  by 
exposing  them  to  a  beam  of  radiant  energy  of  known  intensity  for  a  fixed 
duration.  Time  to  ignition  i.s  reported  as  the  average  over  four  success¬ 
ful  ignitions  and  four  failures. 

Inspection  of  the  arc-image  furnace  test  data  shows  that  the 
ignition  pressure  threshold  level  is  in  the  vicinity  of  one  atmosphere  for 
many  propellants,  as  is  seen  in  Rigures  (2)  and  (3),  But  in  the  case  of 
ANB-3066  LC,  the  level  is  much  lower,  in  the  vicinity  of  4  psia,  (Figures 
(6)  and  (?)).  The  critical  surface  layer  thickness  seems  to  be  proportional 
to  the  average  oxidizer  particle  size.  Normally,  several  oxidizer  grinds 
are  used  in  any  given  propellant  composition,  making  comparison  difficult. 
However,  the  critical  layer  thickness  is  about  3O  microns  for  ANP-2655  AF, 
with  its  relatively  fine  oxidizer  blend,  and  about  100  microns  for  ANB-3066  LC, 
which  has  a  relatively  coarse  oxidizer  blend. 

The  ignitability  data  is  graphed  in  Figures  (2)  through  (7). 

Short  descriptions  of  the  graphs  follow. 
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FiRures  (2)  qnd  (l);  Environmental  pressure  reported  for 


ignition  is  plotted  against  the  time  to  ignition  for  a  variety  of  propel¬ 
lants.  Data  plotted  in  Figure  (2)  was  obtained  at  a  heat  flux  of  70  cal/ 

2  2 
cm  /sec,  while  data  plotted  'a  Figure  (3)  was  obtained  at  20  cal/cr.  /sec, 

'i'he  propellant  surface  preparation  for  these  tests  was  that  of  a  freshly 

cut  surface,  which  has  become  a  standard  for  this  type  of  test, 

Firures  (l)  and  (5);  Data  acquired  on  surface  effects  are  shown 
on  these  two  graphs.  For  each  of  two  types  of  propellant  compositions,  a 
variety  of  test  specimens  were  prepared,  using  different  methods  oi'  surface 
preparation.  The  specimens  were  tested  at  pressure  levels  of  one,  two  and 
three  atmospheres.  Generalities  are  difficult  to  arrive  at  from  the  data 
obtained  However,  surfaces  prepared  by  casting  against  Styrofoam,  and 
then  .stripping  the  Styrofoam  from  the  surface  after  cure,  appear  most 
readily  ignitsbie.  This  may  be  either  a  result  of  a  r.ore  favorable  sur¬ 
face  composition^  or  a  con.sequency  of  an  increase  in  radiation  absorption 
oau.sed  by  extreme  surface  roughness. 

Specimens  having  su"faces  cast  against  RTV-&0,  a  silicone 
rubber  .mold  release  composition.,  were  also  readily  igni  table  compared  with 
the  st.a;idard  cut  surface  preparation.  The  RTV-6O  surface  is  smooth  compared 
with  mat  obtained  from  stripped  Styrofoam.  Figure  (5)  shows  a  remarkable 
f.act  ccncerning  Ai;B-3Cb6  I.C  propellant  composition,  namely  the  near  pressure 
■  nre.as:  fivity  of  its  ign!  tabi  li  ty  over  a  pressure  range  of  one  to  three 
atmospheres, 

FiRures  (13^  and  i~);  The  ignitability  of  one  propellant, 

Aii'B- "’066  LC  having  a  .surface  ra.at  sRainst  RTV-60,  was  found  over  a  wide 


■il- 


'igure  (2 ) 


(ns oc^ 


Figure  ( 


(raee) 


range  of  pressure ^  1.5  psia  up  to  65  psia,  and  at  flux  levels  of  70  and 
40  cal/cEi  /sec.  Specimens  cf  t-he  same  propellant  having  a  fresh  cut  sur¬ 
face  were  also  tested  for  comparative  purposes. 

As  noted  previously;,  this  test  sequence  showed  that  propellant 
having  a  surface  cast  against  RTV-60  is  more  readily  ignitable  than  one 
having  a  surface  prepared  by  cutting.  These  data  also  show  an  ignition 
pressure  threshold  of  about  4  psia  for  this  propellant  composition,  which 
compares  well  with  the  commonly  found  level  of  about  one  atmosphere. 

Figures  (8).  (9^.  and  (10);  Illustrations  of  the  physical 
variations  brought  about  by  various  methods  of  surface  preparation  are 
shown  in  forty  power  photcniicrographs  of  some  of  the  samples  tested. 

Most  specimens,  for  example  those  cast  against  RTV-60,  the  cut  surface, 
and  those  cast  against  155'^F  Polywax,  have  a  smooth  appearing  surface. 
However,  the  surface  cast  against  155°?  Polywax  has  a  distinct  coating  of 
the  wax  remaining  on  the  surface.  Even  though  all  of  these  propellant 
sar-.ples  have  similar  appearances,  the  results  of  the  ignitability  test 
show  a  real  difference.  Samples  of  the  remaining  type  of  surface  prepara¬ 
tion  namely  that  resulting  from  casting  against  styrofoam,  were  jagged 
and  obviously  easy  to  ignite  regardless  of  method  of  foam  removal. 

I  IT,  IGNITER  INDUCED  KOTOR  PRESSURIZATION  AND  GRAIN  IGNITION 

The  combustion  products  from  the  igniter  pressurize  the  port 
volume  of  the  motor  and  transfer  heat  to  the  grain  surface.  In  this 
section  expressions  will  be  found  from  which  the  igniter  induced 
chamber  pressure  may  be  obtained.  And  for  a  particular  case,  that  of  a 
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AJ'IB-3066  ix: 

Cast  ajjainst  HTV, 
strinped  ai'ter  cure 


ANB-3066  T.C 

Polywax, 

stripped  after  cure 


controlled  flov  igniter  vjth  a  propellant  pyrotechnic,  the  expression  for 
heat  fl'ix  to  the  propellant  grain  is  found.  These  tv;o  items,  pressuri¬ 
zation  and  heat  flux,  are  combined  with  ignitability  data  to  provide  an 
analysis  of  2nd  stage  tiinuteman  ignition, 

A  CONSERVATION  EQUATIONS 

Both  mass  and  energy  are  transferi'ed  as  the  combustion  pro¬ 
ducts  of  the  igniter  flow  into  the  port  volume  of  the  motor.  In  this 
section  the.  descriptive  equations  of  mass  and  energy  conservation  will  be 
given.  It  will  be  found  that  the  full  set  of  descriptive  equations  are 
too  cumbersome  to  solve  analytically.  Certain  simplifying  assumptions 
will  be  introduced,  and  the  simplified  equations  solved. 

Mass  Balance:  The  thermodynamic  system  is  the  motor  port 
volume  (combustion  chamber),  and  the  time  interval  is  taken  from  initiation 
of  the  Igniter  up  to  the  instant  of  grain  ignition.  Combustion  of  the 
igniter  py.''otechnic  causes  a  mass  flow  into  the  system.  Simultaneous  with 
pressurization  of  the  port  volume,  there  Is  a  mass  flow  out  the  motor 
nozzle.  The  mtiss  balance  will  be  then,  an  equation  involving  three  terms, 
the  statement  being  as  follows.  The  rate  of  change  of  mass  in  the  system 
equals  the  mass  flow  rate  from  the  igniter  minus  the  mass  flow  rate  out 
the  nozzle.  Symbolically, 

m  =  r'i.  ~  m  (8) 

CIO 

where  m^  =  mass  in  tre  system 

m^  =  mass  from  the  igniter,  assumed  known 

m  =  ma.ss  out  the  nozzle 
o 

The  "dot"  refers  to  tine  differentiation. 
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Flow  out  the  nozzle  will  be  assuned  to  run  full  at  all  tines.  In  other 
words,  the  mass  flow  out  the  nozzle  will  be  given  the  same  form  as  in  the 
sonic  case.  It  can  be  shown  (Reference  (9))  that  when  the  exit  pressure 
is  as  high  as  755K  of  the  chamber  pressure,  the  mass  flow  rate  is  already 
90%  of  maximum  under  sonic  conditions.  Then  using  the  sonic  flow  rate 
form, 


K  P 


m  =  — 
0  ’\f 


VT  /T 
c  o 


where 


K  =  discharge  coefficient,  and 


F  =  igniter  combustion  product  force  constant 
Y  =  isentroplc  constant 
=  motor  throat  area 
=  motor  chamber  pressure 

T^,T^  =  gas  temperature  in  motor  chamber  and  Isochoric 
flame  temperature,  respectively 


The  isochoric  flame  temperature,  T^,  is  of  course,  a  constant. 
Ch, amber  temperature,  T^,  is  expected  to  be  a  variable. 

The  mass  in  the  system  is  obtained  from  the  equation  of  state 
of  a  perfect  gas,  effects  of  co-volume  and  j ncondonsibles  being  accounted 
for  in  the  igniter  combustion  product  force  constant,  F,  (Reference  (10)). 
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Then 
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c  (T  m  )  =  C  T.n.  -  C  T  m 

vat  cc  pii  pco 


(11) 


where  the  additional  terminology  nenns 


C  ,C  =  isochoric  and  isobaric  heat  capacities  of  igniter 


V-  p 


combustion  product,  respectively 


=  flame  temperature  in  the  igniter  combustion  chamber 


If  the  temperature  in  the  igniter  chamber  has  achieved  equili¬ 
brium,  then  Ihe  energy  balance,  Eq.  (11 ),  rewritten  in  terms  of 
familiar  parameters  is, 


where 


dP 


=  C  T 

V  o 


(12) 


Heat  Loss  Rate;  The  heat  loss  rate,  H,  will  be  estimated  on 
the  basis  of  a  convective  heat  flux  from  the  igniter  combustion  product  to 
walls  enclosing  the  combustion  chamber.  In  the  case  of  a  propellant 
igniter  the  convective  type  flux  is  most  significant.  The  Nusselt  relation, 
in  terms  of  igniter  mass  flow  rate  and  a  Prandlt  number  near  unity,  is, 
from  Reference  (11), 


^  =  0.023 


where 


h  =  heat  transmission  coefficient 
d  =  hydraulic  diameter  in  motor  chamber 
p  =  dynamic  viscosity 
k  =  gas  conductivity 


-25- 


m.  =  mass  flow  rate  in  the  chamber,  taken  equal  to  irniter 
mass  flow  rate 

=  flow  area  in  chamber 

The  heat  loss  rate,  H,  heat  flux,  q,  and  heat  transmission 
coefficient  are  related  through, 


H=Aq=Ah(T  -T)  (I3) 

s^  s  c  s 

where  is  the  total  exposed  surface  area  of  the  walls  enclosing  the  com¬ 
bustion  chamber,  and  is  an  average  wall  temperature  during  the  process. 
Average  values  will  ue  used  for  most  thermodynamic  parameters,  which  is 
adequate  since  the  heat  loss  correction  is  a  second  order  effect. 

Several  relationships  will  be  used  to  write  the  heat  trans¬ 
mission  coefficient  in  terms  of  parameters  used  previously.  From  kinetic 
gas  theory  the  gas  conductivity  is  related  to  dynamic  viscosity  and  heat 
capacity  by. 


k 


- 

4 


Then  from  Eq«  (13) 


or 


(9  -  e^) 


(15) 
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where  0^  is  the  ratio  of  average  surface  temperature,  ,  to  the  isochoric 
flame  temperature,  T^,  and. 


r= 


.023 


Substituting  this  expression  into  the  energy  balance. 


(16) 

Eq.  (12), 


results  in, 


V 

P 


m. 


1 


Pm.*®  (e  -  0^)  -  yV^  ka^ 


(17) 


The  mass  balance,  Eq.  (9),  is  repeated. 


V 

P  dt 


(18) 


A  glance  at  the  two  conservation  equations  makes  it  clear  that 
there  are  two  dependent  variables,  igniter  induced  chamber  pressure,  P^, 
and  the  gas  temperature  ratio,  0,  to  be  solved  for  in  terms  of  an  independe 
variable,  time.  Further,  each  of  the  equations  contains  non-linear  terms 
in  a  dependent  variable,  which  in  effect  meMns  the  system  cannot  be  solved 
by  elementary  methods  of  integration. 

Simplifications ;  Since  the  non-linearity  of  the  conservation 
equations  inhibits  obtaining  a  solution  by  elementary  methods,  a  means  of 
linearizing  the  equations  will  be  given  here.  The  method  consists  of 
approximating  the  chamber  temperature,  ,  and  substituting  the  approximate 
value  into  the  energy  balance,  Eq.  (17).  This  will  give  an  approximate 
energy  balance,  linear  in  the  denendent  variable  of  chamber  pressure. 
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Consider  a  total  energy  balance  on  the  igniter  conbustion  pro¬ 


duct, 


C  T.m. 

p  1  1 


C  T  ra. 


p  c  .1 


+  H 


(19) 


where  H  is  the  total  heat  lost  to  the  walls  of  the  enclosure,  and  is 
the  total  igniter  charge  weight.  From  Eq.  (15), 

Pm. (0  -  ej 


Substituting  this  expression  for  total  heat  lost  back  into  Eq.  (19), 
yields  an  approximate  average  value  of  the  temperature  ratio, 


e 

a 


© 

s 


where  6  =  T  /T  ,  constant 

a  c  o' 

The  energy  equations,  Eq.  (17),  reads  then. 


V 

F 


dP 

-TT  =  *  -  yV^  K  a.  P 

dt  e  v  a  to 


which  has  an  integral, 


(20) 


(21) 


(22) 
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where 


lii  =  n,  -  Psi.  6 
e  ilia 

=  ■  f  -ri^*  \ 
b(r)  ) 

=  initial  chamber  pressures 

IV.  APPLICATIONS 

In  the  previous  sections,  collected  arc-image  furnace  ignition  test 
data  has  been  exhibited.  Further,  methods  for  calculeting  the  igniter 
induced  chamber  pressure  and  heat  flux  have  been  found.  However,  there  are 
certain  difficulties  in  obtaining  relations  between  the  arc-image  furnace 
test  data,  the  expression  for  igniter  induced  chamber  pressure,  and  ignition 
requirements  of  an  actual  motor.  The  most  serious  difficulty  is  thet  the 
ignitability  data  is  obtained  at  constent  heat  flux  and  constant  environ¬ 
mental  pressure.  In  an  actual  motor  firing,  the  igniter  induced  chamber 
pressure  is  normally  transient,  and  heat  flux  may  vary  also,  since  it  is 
a  function  of  igniter  mass  flow  rate. 

Most  logically,  the  ignition  test  data  should  be  compared  with  the 
average  igniter  induced  chamber  pressure.  Here,  the  average  is  taken  from 
initiation  up  until  some  time,  t,  rather  than  using  the  average  over  the 
entire  igniter  function  time.  As  far  as  the  heat  flux  is  concernad,  it 
will  be  essentially  constant  for  constant  mass  flow  rate.  And  since  flux 
is  proportional  to  the  four-fifths  power  of  igniter  mass  flow  rate,  i.e., 
a  power  less  than  unity,  the  flux  will  be  more  neerly  constant  than  the 
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mass  flow  rate  Itself.  Hence,  in  many  cases  the  approximation  of  constant 
heat  flux  is  valid. 

The  interpretation  given  to  the  arc-inage  furnace  ignition  test 
data  is  as  fellows. 

A  sequence  of  ignitability  tests  are  conducted.  Over  the  sequence 
of  tests  the  heat  flux  remains  constant.  But  each  test  is  conducted  at 
a  different  level  of  environmental  pressure.  The  result  is  a  curve  in  the 
vs  t^  plane,  where  is  the  pressure  required  for  ignition,  and  t^ 
the  time  to  ignition. 

If  the  curve  P^  vs  t,  where  P^  is  the  time  average  igniter  induced 
chamber  pres.eure,  is  plotted  in  the  same  plane  as  P^  vs  t^,  then  the  inter¬ 
section  of  the  two  curves  should  be  the  time  to  ignition  in  an  actual  firing. 
The  igniter  induced  heat  flux  must  be,  of  course,  equal  to  the  flux  at 
which  the  ignitability  data  was  collected.  These  notions  will  be  illustrated 
by  a  numerical  calculation. 

Time  Average  Igniter  Induced  Chariber  Pressure;  The  average  pressure 
up  to  time  t  is  found  from, 

t 

'  t  /  ’’c'T'l  '^T  («) 

O 

where 

P^  =  time  average  chamber  pressure 
P^  =  chamber  pressure 

The  igniter  induced  chamber  pressure,  P^,  is  given  in  Eq.  (22), 
and  from  this  and  Eq.  (23), 
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b(r)drdV  (24) 


’^ftien  the  igniter  mass  flow  rate  is  constant,  Eq.  (24)  can  be 
evaluated  immediately,  and  is 

^0  -  r  ‘  ^  ^  (1  -  (25) 

Minuteman  Ignition t  The  motor  has  a  grain  configuration  in  the 
shape  of  a  hollow  cylinder  with  fins  in  ihe  forward  section,  as  shown  in 
Figure  (ll).  Ignition  is  by  means  of  a  controlled  flow  igniter,  using  pro¬ 
pellant  os  the  pyrotechnic,  A  total  cf  six  igniter  port  areas  direct  the 
flow  into  the  finned  region,  and  the  areas  of  maximum  heat  flux  will  be 
located  there.  Hence,  time  to  ignition  must  be  calculated  on  the  basis  of 
flow  into  tliat  region.  Average  gas  temperature,  however,  will  be  calculated 
on  the  basis  of  flow  into  the  cylindrical  section,  where  the  flux  is  assumed 

average.  The  hydraulic  diameters  in  the  two  flow  areas  are, 

' 

10  inches .  cylindrical  section 

d  = 

1.5  inches,  finned  section 
Other  motor  parameters  are, 

V  =  13,000  in.^ 

A  =  73  in,"^ 

A.  =  8000  in.^ 
b 

^ ^  -  80  in,^ 
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Figure  (11) 


Igniter  and  pyrotechnic  constants  are, 

=  8  Ibs/sec,  constant 
F  =  3  X  10^  ft-lbs/lb 
Y  =  1.18 

t^  =  0.25  sec,  total  Ainction  tine 

ra^  =  2  11b,  total  charge  weight 
6  2 

p  =  1.5  X  lO”  Ib-sec/ft  ,  dynamic  viscosity 

The  heat  transfer  parameter,  P,  (Eq.  (16)),  evaluating  using  the 
•  ’.ri.r.'iiifclors  listed  above  and  in  the  cylindrical  region,  {d  =  10  inches),  has 
the  value, 

r  =  0.415  (Ibs/sec)'^ 


Average  gas  temperature  ratio,  0  ,  (Eq.  (20),  may  be  obtained  then 

and  is. 


0  =  0.67 

a 


Parameters  occurring  in  the  effective  energy  balance,  Eq.  (21), 
and  their  values  are, 


m  =  in.  -P  A. 6  =  6.53  Ibs/sec 
e  1  '  1  a 

a  =  r  Y  ~\l  s"  '  if  A.  =123  sec  ^ 

V  fl  b 

b  =  ^  =  1810  Ibs/inPsec 

For  initial  pressure  near  zero,  the  time  average  chamber  pressure 


is  from  Eq.  (25) > 


?  =  b  b  (1  _  ^-at) 

c  a  2. 

a  t 


and  using  the  values  given  above, 


T5  T  y  O  120  f-.  “0*123t  \ 

P„  =  14.7  -  -r—  (1  -  e  psia 

c  u 
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0St2  250  msec. 


Heat  flux  for  ignition,  q,  (Eq.  (I3))  is  found  from 


.  r'a  "i”'® 

q  =  - - - 

s 

where  p ’  is  the  heat  transfer  pareuneter  evaluated  in  the  finned 
(d  =  1.5  inches).  For 
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&  0  c  ^  X  on  ^ 


Q  =  1840  Btu/lb 


the  value  of  the  heat  flux  is, 

q  =  257  Btu/ft^  sec -^70  cal/cm^sec 

On  Figure  (12),  an  ignition  schematic  is  plotted  for  the  case  given 
above,  with  the  propellant  grain  consisting  of  ANB-3O66  10.  From  the  arc- 
image  furnace  ignition  test  data,  curves  have  been  obtained  of  pressure 
required  for  ignition  vs  ti-r.e,  for  a  flux  of  70  cal/cra  sec,  for  two  types 
of  surface  preparation  (Figure  (6)).  Recall  in  the  computation  given  above, 
the  he  it  flux  for  ignition  was  found  equal  to  the  ignitability  test  flux. 

By  plotting  the  tert  data  on  the  same  graph  as  the  curve  vs  t  the  time 
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to  ignition  is  found  as  the  intersection  of  the  curves.  Thiii  is  evident 
from  Figure  (12). 

Recapitulating,,  from  motor,  igniter,  and  pyrotechnic  properties, 
the  igniter  ^ndaced  motor  pressure  and  heat  flux  were  found.  The  igniter 
induced  motor  pressure  was  averaged  over  time,  and  compared  with  the  pres¬ 
sure  requirement  for  the  condition  of  propellant  ignition.  The  latter 
curve  is  experimentally  obtained  by  ignitability  measurements  {Figure  (6)), 
at  a  flux  equal  to  that  put  out  by  the  igniter.  When  the  time  averaged 
pressure  curve  equals  the  required  pressure,  ignition  should  occur. 

A  safety  factor  for  motor  ignition  is  defined  as  the  ratio  of  igniter 
function  time  to  ignition  time.  For  the  case  of  2nd  stage  Minuteman,  using 
ANB-3066  LC  propellant  cast  against  RTV,  the  safety  factor  based  on  a 
function  time  of  250  msec  is, 


If  the  surface  is  scraped  or  cut,  the  safety  factor  is  reduced  to  2.3- 
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